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Abstract 

The NPARC Alliance is a partnership between the 
NASA Lewis Research Center (LeRC) and the USAF 
Arnold Engineering Development Center (AEDC) dedi- 
cated to the establishment of a national CFD capability, 
centered on the NPARC Navier-Stokes computer pro- 
gram. The three main tasks of the Alliance are user 
support, code development, and validation. The present 
paper is a status report on the validation effort. It 
describes the validation approach being taken by the 
Alliance. Representative results are presented for lami- 
nar and turbulent flat plate boundary layers, a super- 
sonic axisymmetric jet, and a glancing shock/turbulent 
boundary layer interaction. Cases scheduled to be run 
in the future are also listed. The archive of validation 
cases is described, including information on how to 
access it via the Internet. 

Introduction 

The NPARC Alliance is a partnership between the 
NASA Lewis Research Center (LeRC) and the USAF 
Arnold Engineering Development Center (AEDC) dedi- 
cated to the establishment of a national computational 
fluid dynamics (CFD) capability, centered on the 
NPARC computer program. The NPARC code is based 
on the PARC code, which by 1993 was being widely 
used by a variety of government, industrial, and aca- 
demic institutions. Several of these users had 
approached both LeRC and AEDC about establishing a 
formal organization for the further support, develop- 
ment, and validation of the PARC code. The NPARC 
Alliance was established in 1993 in response to these 
requests. 

The three main tasks of the NPARC Alliance are 
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user support, program development, and validation. 
The Support Team coordinates the release of the soft- 
ware, assists users in its application, resolves bug 
reports and other problems, and manages the various 
communication channels between users and the 
Alliance. The Development Team establishes direc- 
tions for future development of the NPARC code, 
develops enhancements, and incorporates 
improvements contributed by other developers. The 
Validation Team is responsible for validating the 
NPARC code for a wide range of flow parameters and 
geometric configurations, and for establishing an 
archive of cases that can be accessed by the NPARC 
community to support independent assessment of the 
code’s capabilities. 

This paper describes in general the approach being 
taken in the NPARC validation effort. Representative 
results are presented from validation and example cases 
already run, and cases scheduled to be run in the future 
are described. Finally, the validation archive is dis- 
cussed, including information on how to access it via 
the Internet. Additional details on all the cases pre- 
sented here may be found in the documentation avail- 
able in the NPARC validation archive. 

Validation Approach 

The validation effort is intended to establish the 
basis upon which confidence in results produced by 
NPARC is founded, and the practical limits on the accu- 
racy of predictions of flow phenomena pertinent to 
propulsion-oriented flows. Such confidence can only be 
achieved through a continuous process of careful appli- 
cation of the code to a wide range of “unit” and “con- 
figuration-oriented” problems and complete documen- 
tation of results. Here, “unit” refers to problems focus- 
ing on a single phenomenon and simple geometries, 
whereas “configuration-oriented” refers to problems 
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which focus on geometries and flows more representa- 
tive of typical propulsion system components. 

A wide variety of propulsion-oriented validation 
problems representing a mix of flows have been identi- 
fied as candidates for validation cases. Cases are run to 
determine the strengths and weaknesses of the NPARC 
code for a variety of geometric configurations, and over 
a range of flow parameters. Computed results are com- 
pared with benchmark-quality experimental data, well- 
accepted computational results, and/or analytic solu- 
tions. 

Initial validation efforts have focused on “unit” or 
single phenomenon problems, and provide both the ini- 
tial entries into the validation archive as well as serve as 
models for future validation case execution and docu- 
mentation. The following list indicates the type of 
“unit” problems that may be used for NPARC valida- 
tion: 

• Flat plate boundary layers, including heat and mass 
transfer 

• Falkner-Skan flows 

• Curved-wall boundary layers 

• Free shear layers 

• Flow past simple bodies (e.g., a cylinder, sphere, or 
cone) 

• Pressure-driven secondary flows 

• Entrainment flows 

• Shock/boundary layer interactions 

• Separated flows 

• Vortex flows 

• Wake flows 

• Forward and backward facing steps 

The longer term activities are aimed at maintaining 
a process of continuous validation case execution with a 
goal of a minimum of 3-4 validation cases executed per 
year. Increased emphasis is placed on “configuration- 
oriented” problems such as: 

• Propulsive nozzle/afterbody 

• Airbreathing (axisymmetric and 2-D) 

• Rocket (single and multiple nozzles) 

• Airfoil cascade 

• Diffusing duct 

• Propulsion system inlet (subsonic and supersonic) 

• Ejector nozzle 

Model Validation Cases 

The term “validation” has been used in a variety of 
ways in the literature. For this effort, we are guided by 
the following definition, adapted from one given by 
Mehta (1990): 


A code is said to be validated if the following condi- 
tions are met: (1) a comparison of computed results 
with detailed surface and flow field experimental 
data and/or other well-accepted solutions shows 
that the code is able to accurately model the critical 
physics of the flow; (2) the accuracy and limitations 
of the experimental data are known and understood; 
and (3) the accuracy and limitations of the code's 
numerical algorithms , grid density effects, conver- 
gence effects , and physical basis are known and 
understood. The range of applicability of the vali- 
dated code depends on the range of flow parameters 
and/or geometric configurations for which the code 
has been validated. 

Of course, in practice the accuracy and limitations 
of the experimental data and the computational results 
cannot be fully “known and understood ” In addition, 
the degree to which the code must “accurately model 
the critical physics of the flow” will depend on how the 
results are to be used. These factors will inevitably 
introduce some blurring of the line between the states 
of validation and non- validation. Nevertheless, this def- 
inition does serve to provide the necessary philosophy 
that guides the validation effort. NPARC validation 
cases which attempt to meet this strict standard are 
termed model (i.e., ideal) cases. 

Model cases are run to determine the strengths and 
weaknesses of the NPARC code for a variety of geo- 
metric configurations, and over a range of flow parame- 
ters. The Validation Team coordinates with the Devel- 
opment Team to correct weaknesses in the code, and to 
anticipate validation needs relative to on-going and 
planned development work. 

The accuracy and limitations of the code are investi- 
gated by examining the sensitivity of the results to vari- 
ous input options such as mesh density, turbulence 
model, and artificial viscosity model. The validation 
effort is expected to be an on-going activity and, 
through the NPARC Association, users are encouraged 
to propose candidate validation problems and submit 
documentation and results from independent validation 
efforts. 

Each model validation case is documented, in a con- 
sistent format, as part of the validation archive. A con- 
sistent format for the documentation of validation cases 
is necessary to ensure an adequately thorough represen- 
tation of a given case and to permit a comparison of 
conclusions drawn from a variety of cases. The docu- 
mentation highlights the primary focus and pertinent 
findings for each case, and includes sufficient detail to 
allow independent repetition of each case. A Validation 
Abstracts document will be developed summarizing all 
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the model validation cases to compliment the detailed 
information in the validation archive. 

Example Validation Cases 

Example cases are established and documented in 
coordination with the Support team. There are two pri- 
mary goals which the example validation cases are 
designed to meet. The first is to provide users with 
quick, but limited validation of the NPARC software 
over a wide range of flows. These validation cases are 
indicative of the capabilities of the flow simulation pro- 
gram, but do not meet the definition of a model valida- 
tion case in that they do not examine the sensitivity of 
the results to various input options. The second goal of 
the example cases is to provide the new user with clear 
examples of how to properly setup and execute the 
NPARC code for a variety of geometries and flow con- 
ditions. 

The documentation for the example cases includes: 
a description of the problem being solved; a description 
of the computational mesh and initial condition files, 
including, where appropriate, listings of codes and 
input used to create these files; a discussion and listing 
of the input used to run the NPARC code, including 
typical job control commands and the NAMELIST 
input file; and a discussion of the computed results and 
convergence history. The documentation is automati- 
cally provided with the NPARC code as part of the 
NPARC User’s Guide (NPARC Alliance, 1994) and/or 
as a separate document. It is also available as part of 
the validation archive. 

Check Cases 

Check cases will be established to judge the func- 
tionality of a newly installed and/or modified code. 
These will be developed, maintained, and documented 
in conjunction with the Development Team. The pri- 
mary intent of the check cases is to provide the Devel- 
opment team with a tool to ensure the integrity of all 
mechanical aspects of code operation. At least one of 
these cases will be an installation check case that is 
intended for use by new recipients of the NPARC code 
to verify that the code has been properly installed on 
their computer system. 

The documentation for the check cases will be 
included as part of the Developer’s Guide, and will also 
be available as part of the validation archive. 

Flat Plate Boundary Layer 
Laminar Flow 

One of the first model validation cases run with 
NPARC was “incompressible” laminar flow past a flat 


plate. A free stream Mach number of 0.1 was used to 
minimize compressibility effects, and allow comparison 
between the computed results and the exact incom- 
pressible Blasius boundary layer solution. The 
Reynolds number Re x , based on the free stream veloc- 
ity and distance from the leading edge, ranged from 0 to 
200,000. 

Several cases were run for this problem to examine 
the effects of various parameters on the convergence 
rate and the final results. These included the effects of: 
(1) mesh resolution; (2) outer boundary height; (3) 
inflow boundary conditions; (4) inlet plane location; (5) 
initial conditions; (6) artificial viscosity; and (7) using 
the 2D or 3D code. All of the cases were run 12,000 
time steps, with the default value of 2 for IVARDT. 
DTCAP was set equal to 5.0 for the first 3,000 steps, 
10.0 for the next 3,000 steps, and 20.0 for the last 6,000 
steps. For most of the cases, the results were essentially 
the same, and agreed very well with the Blasius bound- 
ary layer solution. Changing the outer boundary height, 
however, did affect the computed values of the normal 
velocity component. 

Three cases were run to investigate the effect of the 
outer boundary height. All used 126 grid points in the x 
direction, evenly distributed between x = -0.25 and 
x = 1.0, where x = 0.0 corresponds to the plate leading 
edge. For the y direction, a Blasius coordinate transfor- 
mation was used, with 

, >. 1/2 

(1) 

where 77 is the Blasius similarity coordinate, u e is the 
free stream velocity and v is the kinematic viscosity. 
First the 77 coordinates were computed, with 21 points 
evenly distributed between tj = 0.0 and 7] - 4.0. Above 
7 = 4.0, which corresponds approximately to the 
boundary layer edge in the Blasius solution, the 77 coor- 
dinate was stretched geometrically using 

V k = TJk-l + r ( 7 lk-l - flk- 2) (2) 

where k is the index in the y direction. For x > 0.25, 
the y coordinates were computed from Equation (1). 
For x < 0.25, the y coordinates were set to those at 
x = 0.25 to prevent the height of the computational 
domain from approaching zero at the plate leading 
edge. 

For the first case, 24 points were added above 
7 ] = 4.0, which gave a maximum k value of 45 and 
7 45 = 23.47. For the second case, 32 points were 
added, which gave 7753 = 48.25. Both of these cases 
used a stretching factor of r = 1.1. For the third case, a 
stretching factor of r= 1.05 was used, and 65 points 
were added, which gave 77 86 = 99.93. (Initially, r = 1.1 
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was also used in the third case, but the solution was 
unstable. It was thought that the large grid spacing in 
the y direction far away from the plate may have been a 
factor in causing the instability, and lowering the 
stretching factor solved the problem.) 

Five boundary segments were used: (1) the 
upstream inflow boundary; (2) the downstream outflow 
boundary; (3) the symmetry plane on the lower bound- 
ary upstream of the plate; (4) the flat plate itself; and (5) 
the upper free stream boundary. The boundary condi- 
tions used, including the code number in the NPARC 
input, are summarized in the following table. 


Boundary 


Bssa 

Upstream 

Fixed conditions 

-10 

Downstream 

Subsonic outflow 

0 

Symmetry plane 

Symmetry 

50 

Flat plate 

No-slip adiabatic 

60 

Upper 

Free stream 

7 


The initial conditions were uniform flow at M = 0.1 for 
all three cases. 

The second-order artificial viscosity coefficient 
DIS2 was kept at its default value of 0.25 for the first 
6,000 time steps, then lowered to 0.0 for the last 6,000 
time steps. The fourth-order coefficient DIS4 was kept 
at its default value of 0.64 for all 12,000 time steps. 

The computed u- and v-velocity profiles in the 
boundary layer region for the three cases are shown in 
Figures 1 and 2. These profiles are at xlL = 0.75, and 
are after 12,000 time steps. Also shown is the exact 
Blasius boundary layer solution (White, 1974). The u- 
velocity results are nearly identical for all three cases. 
However, the v-velocity profiles show significandy bet- 
ter agreement with the Blasius soludon when the height 
of the free stream boundary is increased. (Note, how- 
ever, that in the Blasius similarity soludon the v- 
velocity is normalized differently than the w-velocity, 
and is actually much smaller in magnitude than the w- 
velocity. In the Blasius solution, Re\ n vfu c approaches 
a constant value of 0.8604 far from the plate. At this x 
location on the plate (xlL = 0.75), this corresponds to 
v!u e = 0.00222.) 
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Figure 1. Effect of free stream height on u-velocity 
profiles for laminar flow past a flat plate. 



Figure 2. Effect of free stream height on v-velocity 
profiles for laminar flow past a flat plate. 

Itirbulent Flow 

Another basic model validation case run with 
NPARC was turbulent flow past a flat plate. A free 
stream Mach number of 0.2 was used to minimize com- 
pressibility effects, and allow comparison between the 
computed results and incompressible experimental data. 
The Reynolds number Re x , based on the free stream 
velocity and distance from the leading edge, ranged 
from 0 to 1 x 10 7 . 
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As in the laminar flow problem, several cases were 
run to examine the effects of various parameters on the 
convergence rate and the final results. These included 
the effects of: (1) mesh resolution; (2) the y + value at 
the first grid point away from the wall; (3) turbulence 
model; (4) inflow boundary conditions; (5) inlet plane 
location; (6) artificial viscosity; and (7) the time step 
selection option. In general, the computed results 
agreed very well with experimental data. The most sig- 
nificant differences between calculations were due to 
inadequate mesh resolution, too-large y + values, and the 
choice of turbulence model. As an example, the follow- 
ing discussion presents results for different y + values, 
with the Baldwin-Lomax (1978) turbulence model. 


The computational mesh was generated alge- 
braically. Grid points in the x direction were packed 
near the leading edge of the plate, and in the y direction 
near the plate surface, using a Roberts transformation 
(Roberts, 1971). 

For the x, or j, direction, 126 points were used, with 
the leading edge at x = 0.0 corresponding to j = 26. 
The computational domain extended from 
X] =-1.0352, upstream of the leading edge, to 
X 12 6 = 10.0. 

For the y, or k, direction, 76 points were used, with 
51 in the boundary layer. The boundary layer thickness 
was estimated as follows (Daily and Harleman, 1966): 


5 = 


0. 38x 
Re x ' 15 


(3) 


Different amounts of grid packing were used within the 
boundary layer to give different y + values at the first 
grid point away from the wall. For y > 5, the y coordi- 
nates were stretched geometrically using 


y* = y*-i + Ky*-i - y*- 2) ( 4 ) 


The stretching factor r was 1.1. For j < 36, the y coor- 
dinates were set to those at j = 36 to prevent the height 
of the computational domain from approaching zero at 
the plate leading edge. 

It should be noted that since the number of points in 
the boundary layer was fixed, tighter spacing near the 
wall meant larger spacing near y = S. Since there were 
25 points outside the boundary layer for all the cases, 
the height of the outer boundary varied with the degree 
of packing near the wall. Also, since y + depends on the 
computed flow field in addition to the physical y coordi- 
nate, the y + values did not scale exactly with the y 2 val- 
ues. The values of some of these grid-related parame- 
ters are shown in the following table. 


Case 

yi/s 

y + 

3Wcr ^ 

8 

0.82 XlO -4 

0 . 0 - 0.5 

10.070 

1 

1.64 xlO -4 

0 . 0 - 1 .0 

9.216 

5 

3.33 X 10 -4 

0.0-1 .9 

8.329 

6 

6.56 x 10 -4 

0.0-3 .4 

7.468 

7 

13.11 xlO -4 

0.0-5.5 

6.578 

24 

26.29 xlO -4 

0.0-8 .2 

5.678 


As in the laminar case, five boundary segments 
were used: (1) the upstream inflow boundary; (2) the 
downstream outflow boundary; (3) the symmetry plane 
on the lower boundary upstream of the plate; (4) the flat 
plate itself; and (5) the upper free stream boundary. 
The boundary conditions were the same as in the lami- 
nar case. The initial conditions were uniform flow at 
M = 0.2. 

Cases 1, 5-7, and 24 were run 14,000 time steps 
with DTCAP = 5.0. Because of the tighter mesh spac- 
ing, however, case 8 required more iterations and a 
smaller time step for stability. It was first run 2,000 
steps with DTCAP = 5.0 to get past the starting tran- 
sient. During this initial run, the parameter PCQMAX, 
which was left at its default value of 10.0, controlled 
the actual time step size. The value of DTCAP was 
then lowered to 1.0 for 2,000 iterations, then raised to 
2.0 for the rest of the computation. A total of 28,000 
iterations were taken for case 8. 

The computed w- velocity profiles for cases 1, 6-8, 
and 24 are shown in Figure 3, along with the experi- 
mental data of Weighardt (Coles, D. E., and Hirst, E. 
A., 1968). The results for the three cases with y + < 3.4 
are essentially the same, but differences can be seen for 
cases with larger y + values. As y + increases, the veloc- 
ity profiles tend to become less steep. 



Figure 3. Effect of initial y + on n-velocity profiles for 
turbulent flow past a flat plate, Baldwin-Lomax model. 

The computed skin friction coefficients are com- 
pared with the experimental data of Wieghardt, and 
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with values from two different correlation formulas 
(White, 1974), in Figure 4. The results for the two 
cases with y + everywhere below 1.0 are essentially 
identical. As y + increases, the skin friction decreases, 
which is consistent with the velocity profile results 
shown in Figure 3. 
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o 0.006 


•c 0.004 


0.002 


0.0x1 0 7 0.2x1 0 7 0.4x1 0 7 0.6x10 7 0.8x1 0 7 I.OxlO 7 

Distance, Re x 

Figure 4. Effect of initial y + on skin friction distribu- 
tion for turbulent flow past a flat plate, Baldwin-Lomax 
model. 




NPARC, y* = 0.0-0.5 

— 

NPARC, y + = 0.0-1 .0 

— 

NPARC, y + = 0.0-3.4 

— 

NPARC, y + = 0.0-5.5 


NPARC, y+ = 0.0-8.2 

0 

Experiment, Wieghardt 

□ 

Correlation, c ( = 0.0592Re x * 1/s 

A 

Correlation, C| = 0.455 / ln 2 0.06Re x 




Supersonic Axisymmetric Jet 

Another model validation case currently being com- 
puted is the turbulent supersonic axisymmetric jet flow 
studied experimentally by Eggers (1966). One of the 
runs was used as an example validation case, and some 
of the results from that calculation are presented here. 
This case is described in detail in the NPARC User’s 
Guide (NPARC Alliance, 1994). 

The jet is produced by an axisymmetric convergent- 
divergent nozzle. This study focuses primarily on the 
development of the jet downstream of the nozzle exit, 
although care is taken to properly model the flow within 
the nozzle. The basic geometric configuration is shown 
in Figure 5. 


Jet Mixing 



Centerline 


Figure 5. Geometric configuration for axisymmetric jet 
flow. 

The nozzle total pressure and temperature were 
162.2 psia and 525 °R, respectively, and these values 


were selected as the reference values. The reference 
length X r was 1 ft, giving a reference Reynolds number 
Re r = a r X r lv r of 78.8 x 10 6 , where a r and v r are the 
speed of sound and kinematic viscosity evaluated at the 
reference conditions. The nozzle was operated at a 
pressure ratio of 0.09063, corresponding to perfect 
expansion with an exit Mach number of 2.22. The noz- 
zle exit radius was 0.5035 in. 

The computational mesh was generated using 
GRIDGEN (Steinbrenner, J. P., Chawner, J. R., and 
Fouts, C. R., 1990). The mesh consisted of 336 points 
in the x, or y, direction, and 107 points in the r, or k, 
direction. In the nozzle interior, there were 188 and 59 
points in the x and r directions respectively. The com- 
putational domain extended approximately 145 nozzle 
exit radii downstream of the nozzle exit. In the radial 
direction, the grid extended approximately 12 nozzle 
exit radii at the nozzle exit plane, and approximately 48 
nozzle exit radii at the exit of the computational 
domain. The grid in the region of the nozzle is shown 
in Figure 6. For clarity, only every fourth grid line in 
the x direction and every sixth grid line in the r direc- 
tion is shown. 



flow. 

As mentioned earlier, care was taken in constructing 
the grid interior to the nozzle to ensure adequate repre- 
sentation of the nozzle exit flow. It was found through 
preliminary inviscid studies that the quality of flow at 
the exit of the nozzle was highly dependent on the axial 
spacing of the grid in the diverging section. In particu- 
lar, adequate nozzle wall resolution was required to 
minimize the generation of a series of compression 
waves from the discrete representation of the wall 
geometry. For this problem, 110 points in the x direc- 
tion in the diverging section were found to be adequate. 

As an additional note, since the original report by 
Eggers provides no details as to how the nozzle was 
connected to the supply reservoir, the inviscid study 
also focused on nozzle entrance effects. The nozzle 
was assumed to draw from both a constant area duct as 
well as an “infinite” radius volume. This study 
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indicated no effect of the reservoir geometry on nozzle 
wall pressure distribution downstream of the throat. 
Thus, for the purpose of this example case, the nozzle 
was assumed attached to a constant area duct. 

Since the focus of this example is the development 
of the mixing region downstream of the nozzle exit, 
grid resolution within the jet shear layer throughout its 
development region was a principal concern. The mini- 
mum radial grid packing is driven by the shear layer 
thickness at the nozzle exit. To estimate this thickness 
and also provide guidance relative to the radial grid 
spacing adjacent to the nozzle wall, the finite difference 
boundary layer code BLAYER (Hodge, B. K., and 
Adams, J. C., Jr., 1978) was used to estimate the bound- 
ary layer growth within the nozzle. This required as 
input a nozzle wall pressure distribution, which was 
taken from the nozzle entrance effects study. From 
these results, the physical grid spacing required to 
maintain at least one grid point in the laminar sublayer 
(i.e., y + < 1) was estimated to be 1 x 1(T 5 inches. The 
^-constant grid line emanating from the nozzle exit lip 
(i.e., k = 59) was positioned to approximate the center 
of the jet shear layer based upon Eggers’ data, and pro- 
vided the location about which radial grid packing for 
shear layer resolution was centered. This radial pack- 
ing was gradually relaxed toward the flow domain exit 
boundary where the radial grid spacing was uniform 
between the centerline and the k = 59 grid line and 
gradually increased from k = 59 to the outer boundary 
at k = 107. The clustering function used was the 
default hyperbolic tangent function in GRIDGEN. 
Points were also packed in the x direction near the noz- 
zle exit lip, where the first axial grid spacing down- 
stream of the lip was 8 x 10” 3 inches. 

A total of eight boundary segments were used in the 
NPARC calculation: (1) the nozzle inflow boundary; (2) 
the downstream outflow boundary; (3) the nozzle exit 
lip; (4) the upstream free field inflow boundary (exterior 
to the nozzle); (5) the centerline; (6) the upper 
ffeestream inflow boundary; (7) the nozzle interior wall; 
and (8) the nozzle exterior wall. The boundary condi- 
tions used, including the code number in the NPARC 
input, are summarized in the following table . 3 


Boundary 

Condition 

Code# 

Nozzle inflow 

Subsonic inflow 

0 

Downstream outflow 

Subsonic outflow 

0 

Nozzle exit lip 

Slip wall 

50 

Upstream free field 

Subsonic inflow 

0 

inflow 



Centerline 

Axis of symmetry 

51 

Upper ffeestream 

Subsonic inflow 

0 

inflow 



Nozzle interior wall 

No-slip adiabatic 

60 

Nozzle exterior wall 

Slip wall 

50 


The initial conditions were established by setting 
the temperature and pressure equal to the ambient val- 
ues, and each velocity component to zero for all points 
in the computational domain except those interior to the 
nozzle. There, total conditions were set to the supply 
values, and velocities were initialized by linearly inter- 
polating in the x direction using the j index, from 
M = 0.3 at the nozzle inflow plane to a M = 2.22 at the 
nozzle exit plane. 

A total of 10,000 time steps were taken, with the 
default value of 2 for IVARDT. DTCAP was set to 5.0 
for the first 1500 steps, and 0.4 thereafter. For the 
entire solution, the maximum change allowed in pres- 
sure or density over a time step was 25%. The artificial 
viscosity coefficients DIS2 and DIS4 were left at their 
default values of 0.25 and 0.64 for the first 1500 steps. 
DIS2 was lowered to 0.12 for the next 3500 steps, and 
0.0 thereafter. The problem was run inviscidly for the 
first 1000 steps to quickly develop the initial plume 
characteristics. A combination Baldwin-Lomax (1978) 
and Thomas (1979) algebraic turbulence model was 
used for the next 1 100 steps to initialize the turbulent 
quantities. The Chien k-e model (1982) was used for 
the rest of the calculation. 

Figure 7 presents the computed nozzle centerline 
velocity distribution compared with the measurements 
of Eggers (1966). The local jet centerline velocity is 
nondimensionalized by the nozzle exit velocity and dis- 
tances are measured relative to the nozzle exit station 
and are scaled by the nozzle exit radius. 


3. Daring problem start-up, the nozzle interior wall was actually 
specified as a slip wall (i.e., code #50), then changed to a no- 
slip adiabatic wall after 1000 time steps. 





Figure 7. Nozzle centerline velocity. 

The computed centerline velocity agrees well with the 
experimental data. As can be seen the figure, closure of 
the jet core occurs at approximately 20 radii. For refer- 
ence, from the calculations the nozzle centerline veloc- 
ity first reaches a subsonic value at approximately 45 
radii and flow conditions become “incompressible” 
(i.e., the centerline Mach number drops below 0.3) at 
approximately 120 radii. 

In Figures 8 and 9, the computed velocity profiles 
are compared with experimental data at two locations 
downstream of the nozzle exit. The experimental data 
were taken at x!r e = 73.80 and 121.30. These positions 
are indicated by the filled circles in Figure 7. The com- 
putational results are at x/r e = 73.67 and 121.46, which 
correspond to the j indices closest to the experimental 
data locations. 




Local u Velocity / Centerline u Velocity 


Figure 9. Velocity profile at x/r e = 121.3. 

The velocity profile at xlr e = 73.80 agrees well with the 
experimental data, but the profile at x/r e = 121.3 does 
not agree as well. This may be due to the influence of 
the extrapolation boundary condition used at the nearby 
exit plane on the computed results in this region. In 
addition, the x-momentum flux for this case was contin- 
uing to decrease when the calculation was stopped 
(NPARC Alliance, 1994), and continuing the calcula- 
tion may improve the agreement. 

Glancing Shock/Boundary Layer Interaction 

Another model validation study currently underway 
is the interaction of a glancing shock wave with a turbu- 
lent boundary layer on a flat plate. One of the cases is 
being used as an example validation case, and some of 
the results from that calculation are presented here. 

The glancing shock/turbulent boundary layer inter- 
action is produced by a Mach 3.0 flow past a sharp 10° 
wedge, or fin, mounted on a flat plate. Figure 10 shows 
the geometric configuration and the outline of the com- 
putational domain. The location of the inviscid shock is 
also shown in the figure, with a shock angle of 29.5°. 
This flow has been studied experimentally by Kim, et. 
al. (1990, 1991), and tabulated experimental values of 
skin friction coefficient, surface flow angle, and surface 
static pressure have been published by Settles and Dod- 
son (1991). The Reynolds number in the experiment, 
based on free stream conditions, was 6.19 x 10 7 /m, and 
the free stream total pressure and temperature were 827 
kPa and 294 K, respectively. 
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Figure 10. Geometric configuration for the glancing 
shock/turbulent boundary layer interaction. 

Cartesian coordinates are defined with the origin at 
the fin leading edge, and the flat plate lies in the x-z 
plane. The measured boundary layer thickness 38 mm 
upstream of the fin leading edge was 3.02 mm, and this 
is used as the reference length. Lengths shown in the 
figure are non-dimensional. The leading edge of the flat 
plate was at x = -71.523. 

In the experiments, data were taken along a circular 
arc of radius R, centered at the fin leading edge. For the 
skin friction coefficient and surface flow angle 
R = 37.85, and for the surface static pressure 
R = 33.64. The angular location along this arc is 
denoted by p, as shown in Figure 10. 

The 3-D computational mesh was generated alge- 
braically. In the x direction, 51 points were used, 
evenly spaced. In the y direction, normal to the flat 
plate, the first interior point was placed at the location 
corresponding to y + = 0.5, where y + was computed 
assuming the skin friction coefficient cy = 0.00152. 4 
Points were then added, stretched geometrically using 

y k = yk- 1 + r y (y k -\ - y k - 2) ( 5 ) 

The stretching factor r y was 1.2. Equation (5) was used 
until 

A y = y k - y k -\ > (A y)max ( 6 ) 

where (Ay)^ = 2000. Above this point an evenly 
spaced mesh was used. The same y-mesh was used for 
all x and z. A similar procedure was used in the z direc- 
tion, intersecting the fin. At the experimental measure- 
ment station for static pressure, R = 33.64, the first inte- 
rior point was placed at the location corresponding to 
z + = 0.5, where z + was computed assuming 
Cf = 0.002. 5 Points were then added, stretched geomet- 
rically by the factor r z = 1.5, until the spacing 
A Z > (A z)max> where (A z)max = 5000. Beyond this point 

4. This is the experimental value at the incoming boundary layer 
measurement station, 38 mm upstream of the fin leading edge. 

5. This is a rough estimate of the value on the flat plate, near the 

fin, based on the experimental data along the data measurement 
radius shown in Figure 10. 


an evenly spaced mesh was used. This z-distribution 
was then repeated, proportionally, for all x and y. Near 
the x = 0 station the outer z boundary, which would oth- 
erwise have had a slope discontinuity corresponding to 
the fin leading edge, was smoothed using a fifth-order 
polynomial. The sharp fin leading edge at the inner z 
boundary, of course, was retained. 

Seven boundary segments were used: (1) the 
upstream inflow boundary; (2) the downstream outflow 
boundary; (3) the flat plate surface; (4) the upper free 
stream boundary; (5) the symmetry plane at z = 0 
upstream of the wedge; (6) the wedge surface; and (7) 
the outer z boundary. The boundary conditions used are 
summarized in the following table. 


Boundary 

Condition 

Code# 

Upstream 

Fixed conditions 

-10 

Downstream 

Extrapolation 

3 

Flat plate 

No-slip adiabatic 

60 

Upper 

Extrapolation 

3 

Symmetry plane 

Symmetry 

50 

Wedge 

No-slip adiabatic 

60 

Outer 

Extrapolation 

3 


The conditions at the upstream inflow boundary 
were computed by interpolation from a 2-D NPARC 
calculation of the turbulent boundary layer flow past the 
flat plate upstream of the wedge. The initial flowfield 
was set equal to the flow at this upstream boundary. 

The default Baldwin-Lomax algebraic turbulence 
model was used for the first set of runs. Additional runs 
were then made with the Chien k-e model, using the 
converged Baldwin-Lomax results as initial conditions. 
For the Baldwin-Lomax calculation, DTCAP was set 
equal to 5.0 for the first 800 time steps, 1.0 for the next 
800 steps, and 0.5 for the rest of the calculation. For 
DTCAP > 1.0, the actual time step was computed by 
the code to limit the maximum change in density or 
pressure to 10%, the value specified by PCQMAX. A 
total of 14,400 time steps were taken. The Chien k-e 
calculation also used a DTCAP value of 0.5, for 7500 
additional time steps. 

Figure 1 1 shows computational “oil-flow” patterns 
on the flat plate. 6 Starting at the wedge leading edge 
station, “particles” were released at every grid point 
along the outer z boundary, selected grid points along a 
line at x = 0, and at every grid point along a line near 
the inner z boundary. 


6. Computationally, these are actually in the plane of the first grid 
point above the surface, since all the velocities on the surface 
are zero. 
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Figure 11. Computational “oil-flow” patterns, glanc- 
ing shock/turbulent boundary layer interaction. 

The computed values of the skin friction coefficient, 
the surface streamline angle, and the static pressure are 
compared with experimental data reported by Settles 
and Dodson (1991) in Figures 12(a) through 12(c). 
Computed results are shown for both the Baldwin- 
Lomax and Chien k-e calculations. The skin friction 
data were obtained using a laser interferometer skin 
friction meter, which measures the thickness of an oil 
film on the surface (Kim and Settles, 1990). The error 
bars in the skin friction plot are those reported by Set- 
tles and Dodson (1991), and reflect the repeatability of 
the measurement. The experimental surface streamline 
angles were measured from kerosene-lampblack flow 
visualization patterns, with an accuracy of ±5%. The 
static pressure data were measured using taps on the 
surface, and are believed accurate to within ±3%. Data 
points without error bars (and this includes all the sur- 
face streamline angle and pressure data) have errors 
comparable to, or less than, the size of the symbols in 
the plot. 





Figure 12. Computed results compared with experi- 
mental data, glancing shock/turbulent boundary layer 
interaction. 


Additional Validation Cases 

Work is currently planned or underway on a variety 
of additional model and example validation cases. 
Model validation studies currently scheduled to be run 
over the next two years include: 

• Turbulent flow past a flat plate with and without 
heat transfer 

• Subsonic diffuser flow (Dudek, Georgiadis, and 
Yoder, 1996) 

• Subsonic flow past an airfoil 

• Glancing shock wave/turbulent boundary layer 
interaction 

• Rearward facing step 

• Subsonic flow through an S-duct 

• Supersonic free jet flow 
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Selected cases from some of the above studies will also 
be documented as example validation cases. Additional 
example cases currently scheduled include: 

• External flow past an aft boat-tail configuration 

• Unsteady normal shock, in conjunction with a 
NASA LeRC experiment 

• Transition duct 

• Counter-flowing jet 

• Crossed glancing shocks with and without bleed, in 
conjunction with a NASA LeRC experiment 

Data Archive 

For validation cases to be of maximum benefit to 
the user community, the results must be made readily 
available. The information should be also be detailed 
enough to permit the calculations to be repeated with 
relative ease by independent code users. A formal elec- 
tronic archive system has been established to meet these 
criteria. 

This archive is the central repository for the model 
validation cases, example validation cases, and check 
cases computed by and for the NPARC Alliance. It is 
intended for use by NPARC users and developers, and 
allows easy access to the results of latest validation 
studies. The archive includes all the input files used to 
run the cases, the output files, the experimental data 
used for comparison with computed results, and written 
documentation providing an overview of each case and 
a discussion of the results. The archive is accessible 
over the Internet via anonymous ftp at 
info.amold.af.mil/pub/nparc, and from the existing 
NPARC WWW home page located at 
http://info.amold.af.mil/nparc/index.html. 
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